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Abstract 

Lifetime values for alpha decay in even-even nuclei with Z = 84 — 98 and 
N = 128 — 152 have been calculated in the superasymmetric fission model. 
The interaction between the alpha particle and the daughter nucleus has been 
formed in the double folding approach using a density dependent NN inter- 
action. The densities have been obtained using the Relativistic Mean Field 
formalism. The spectroscopic factors for the decays have been deduced and 
are shown to vary smoothly as a function of effective numbers of valence nu- 
cleons, N p and N n chosen with a suitable core. The implication of such a 
smooth behaviour has been discussed. 

PACS : 21.60.Jz, 23.60.+e 

Keywords: Relativistic Mean Field, Superasymmetric fission model, alpha decay 
spectroscopic factor, N p N n scheme 

Simplified parametrization of various nuclear quantities may be obtained if the 
quantities are plotted as a function of N p N n , the product of effective number of 
valance particles (or holes) [lj. Various quantities such as deformation and B(E2) 
valuespj [31 Hj, rotational moments of inertia in low spin states in the rare earth 
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region [5], ground band energy systematics [5], core cluster decomposition in the 
rare earth region[7| and properties of excited states [HI E] have been found to follow 
certain simple trends when expressed as a function of simple product of N p and 
N n or certain simple functions of the above two numbers. Essentially this simple 
functions are seen to represent n — p interaction and bear smooth relationships with 
the observables. 

It has often been pointed out that in view of the change in magic number and shell 
structure in various mass regions, the conventional counting of valence protons and 
neutrons may be inadequate. Extraction of effective number of valence particles in 
the N p N n scheme may significantly improve the predictive capability of the scheme as 
well as point to the emergence of new shell structure in different mass regions[Tl[TU]. 
For example, Zhao et a/.[TU] have used the concept of effective valence number of 
protons and neutrons to study the breakdown of the shell gap at Z=64 in even-even, 
odd and odd-odd nuclei. Wolf and CastenJTTJ also have studied the shell gap at 
Z = 64 for N p N n scheme. 

Alpha decay in the heavy actinide nuclei presents a unique opportunity of testing 
various nuclear structure theories. It is known to take place through tunnelling of 
the potential barrier by the alpha particle. The spectroscopic factor in a-decay 
was introduced to incorporate the preformation probability. It contains the nuclear 
structure effects, and may be thought as the overlap between the actual ground 
state configuration of the parent and the configuration described by one a-particle 
coupled to the ground state of the daughter. Thus the tunnelling probability of the 
alpha particle is a quantity which is expected to be crucially dependent on the n — p 
interaction and the quantity spectroscopic factor may show a certain correlation 
when expressed as a function of N p and N n . 

In the present work, we have followed the microscopic superasymmetric fission 
model. The interaction potential between the alpha particle and the daughter nu- 
cleus has been obtained in the double folding model by folding the densities of the 
alpha particle and the daughter nucleus using some suitable interaction. Usually 
the densities are obtained from phenomenological description. However, microscopic 
densities obtained from mean field approaches may be expected to provide a better 
description of the densities and hence that of the process of alpha decay. 

Relativistic Mean Field (RMF) approach is now a standard tool in low energy 
nuclear structure. It has been able to explain different features of stable and exotic 
nuclei like ground state binding energy, deformation, radius, excited states, spin- 
orbit splitting, neutron halo, etc. For a detailed discussion, readers are referred 
to the refs. [T2], [131 El- Later related developments in studies of ground state 
structure include incorporating the effect of the resonant continuum in drip line 
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nuclei in both mean field [TB] and Hartree Bogoliubov approach |lb'j. study of new 
Lagrangian densities with coupling between mesons [TT]. etc. It is well known that in 
nuclei far away from the stability valley, the single particle level structure undergoes 
certain changes in which the spin-orbit splitting plays an important role. Being 
based on the Dirac Lagrangian density, RMF is particularly suited to investigate 
these nuclei because it naturally incorporates the spin degrees of freedom. 

There exist different variations of the Lagrangian density as well as a number of 
different parameterizations in RMF. Recently, a new Lagrangian density has been 
proposed [IT] which involves self-coupling of the vector-isoscalar meson as well as 
coupling between the vector-isoscalar meson and the vector-isovector meson. The 
corresponding parameter set is called FSU GoldfTT]. This Lagrangian density was 
earlier employed to obtain the proton nucleus interaction to successfully calculate 
the half life for proton radioactivity [18]. We have also applied it to study alpha 
radioactivity in superheavy nuclei with A > 282 [19] and cluster radioactivity[20] in 
the same approach followed in the present work. This density seems very appropriate 
for a large mass region viz. medium mass to superheavy nuclei. In this work also, 
we have employed FSU Gold. 

In the conventional RMF+BCS approach for even-even nuclei, the Euler-Lagrange 
equations are solved under the assumptions of classical meson fields, time reversal 
symmetry, no-sea contribution, etc. Pairing is introduced under the BCS approxi- 
mation. Since accuracy of the nuclear density is very important in our calculation, 
we have solved the equations in co-ordinate space. The strength of the zero range 
pairing force is taken as 300 MeV-fm for both protons and neutrons. 

The microscopic density dependent M3Y interaction (DDM3Y) is obtained from 
a finite range nucleon interaction by introducing a density dependent factor. This 
class of interactions has been employed widely in the study of nucleon-nucleus as 
well as nucleus-nucleus scattering, calculation of proton radioactivity, etc. In this 
work, we have employed the exponential density dependent interaction DDM3Y1 

v(r, Pl , p 2 , E) = C(l + a exp {-(3{ Pl + p 2 ))(l - 0.002E)u M3Y (r) (1) 

used in Ref. [21] to study alpha-nucleus scattering. It uses the direct M3Y potential 
u M3Y (r) based on the G-matrix elements of the Reid[22] NN potential. The weak 
energy dependence was introduced [23j to reproduce the empirical energy dependence 
of the optical potential. The parameters used are the standard values viz. C = 
0.2845, a = 3.6391 and (3 = 2.9605fm 2 . Here pi and pi are the densities of the a- 
particle and the daughter nucleus and E is the energy per nucleon of the a-particle in 
MeV. This interaction has been double folded with the theoretical densities of alpha 
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particle and the daughter nucleus in their ground states using the code DFP0T[24J. 
The assault frequency has been calculated from the decay energy following Gambhir 
et a/.[25]. 

We have studied the alpha decays between the ground states of even-even nuclei 
for parents with Z = 84 — 98 and N = 128 — 152. The spectroscopic factor has been 
calculated as the ratio of the calculated half life to the experimentally observed value. 
Obviously it is expected to be less than unity. The results for the spectroscopic 
factors are presented in Table [TJ Theoretical calculations exist for spectroscopic 
factors for some of these nuclei. Our results are comparable to those values. For 
example, our calculated value for S of 212 Po is 0.019 as compared to theoretical value 
0.025 deduced in [22]. A value of 0.031 was obtained by Mohr[27] in a double folding 
model calculation using density from experimentally known charge distribution. 

The spectroscopic factors for different chains of isotopes generally follow a trend 
though there are exceptions. However, we find that expressed in terms of the Casten 
factor P = N p N n / (N p + N n ), the points, with some exceptions, show a smooth trend 
when the N p and N n values are, or in other words, the core is, chosen suitably. As 
earlier pointed out, this method has already been followed [Tj, QUI- However, the 
implications of such choices of N p and N n values are usually not always obvious. 
We will invoke a number different of subshells in various nuclei that are designed 
to produce a smooth trend and then give microscopic justifications for these choices 
later in the paper. 

For Po(Z = 84) and Rn(Z=86) isotopes, the masses are very close to the doubly 
magic core Z = 82, N = 126. Thus N p and N n are chosen simply as N p = Z — 82 
and N n = N — 126, respectively. However, the values for Ra(Z = 88) isotopes fall on 
the smooth curve only if we choose N p = 4. This can be justified on the assumption 
of Z = 92 as a semiclosed shell as explained later. Similarly, for Th(Z = 90) nuclei, 
the effective value needed to be chosen is also N p = 4. U(Z = 92) nuclei present a 
somewhat different scenario. For the two lightest U nuclei studied, i.e. 224 > 226 U, we 
choose N p = 2 while for all the heavier isotopes, N p — 4 is taken as the effective 
valence proton number. The valence neutron number is taken to be N n = N — 126 
in all the above nuclei. For Cm(Z = 96) and Cf(Z = 98) isotopes, we choose N p = 2 
and 4, respectively. For Cm isotopes, the valence neutron number is taken to be 
N n = N — 138 while For Cf isotopes, it is taken as N n = 152 — N. However Pu 
(Z = 94) nuclei could not be fitted in the scheme. Interestingly, this is analogous 
to the case of some nuclei with N=90-92 which could not be fitted to the N p N n 
scheme for B(E2) values [TO] and for the ratio between the excitation energies of the 
first 4 + and 2 + states [8]. The resulting plot for spectroscopic factors as a function 
of P(> 1.0) is shown in Figured! The smooth line is drawn only for the purpose of 
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guiding the eye. The results for Pu have also been plotted in Figured! As results for 
spectroscopic factor of Pu did not fall on the smooth line in Figure [T], it is difficult 
to assign N p and N n values in the present procedure. We have assigned N p =2 and 
N n = N — 126 for all the Pu nuclei. 

It will be interesting to trace the origin of the effective number of valence par- 
ticles in different isotopes. As already explained, for Po and Ra nuclei, the number 
of valence particles is simply calculated as those outside the g2 8 Pb core. Our cal- 
culation suggests that the next proton single particle level is 71T/19/2. This fact 
is also experimentally verified by the ground state spin-parity of the odd proton 
nuclei with Z = 83 — 91. Particularly, the ground states of the odd mass nu- 
clei 187 - 215 Bi(Z = 83), 197 - 217 At(Z = 85), 20i-2i 9 p r(z = 207-219 Ac(z = g9)) 
2i3-22ip a ^ = gi) a]} have spin-parity 9/2~, whenever known. Theoretically, we 
find an interesting situation. The level just below the 82 shell gap, i.e. tt3si/2 and 
the three levels above it, i.e. irlhg/2, ^2/7/2 and 7rlzi3/2 are shown in Figure [2] for 
nuclei with N — Z = 44. It can be easily seen that the nlhg/2 state is going down 
sharply with increasing Z, thus opening up another gap at Z=92. Earlier mean field 
calculations of Rutz et al. [2B] also suggested a shell gap at Z = 92 in the vicinity of 
nuclei with N = 126. Hence, in Ra nuclei protons are hole-like and N p is taken as 
4. 

In nuclei with higher Z, there is a possibility of another gap opening up near 
Z = 94. Nuclei with Z = 84 — 88 that we have studied have small deformation in 
their ground state. This is evident from the measured B(E2) values. However, near 
Z = 92, nuclei have large deformation. There is a possibility of new gaps opening 
up in this region for deformed nuclei. We should mention here that deformed shell 
gaps are not usually employed in the N p N n scheme. In an earlier calculation with 
RMF, Long et al. (29] have shown that in nuclei with deformation (3 ~ 0.3 and with 
Z and iV near 96 and 144, respectively, gaps of the order of 1.5-2 MeV appear at 
Z = 94, 96. There is a similar example in the neutron sector. Rare-earth nuclei such 
as Hf, Lu, etc are known to exhibit a deformed shell gap at N = 94. This possibly is 
the reason that Cm and Cf isotopes have N p = 2 and 4, respectively while Th with 
Z = 90 has 4 proton holes. However, our calculation does not include the effects of 
deformation and it is difficult to be absolutely certain of this explanation. This fact 
also also may explain why 224 > 226 U have N p = 2 but fails for the heavier U nuclei 
where it was necessary to assign N p — 4. 

A similar situation also exists in the neutron sector. In nuclei with neutron 
number N > 126, the first level beyond N = 126 is calculated to be the uliu/2 
state. In lighter nuclei this state is close to the other states in the same shell. 
However, as neutron number increases, it comes down in energy rapidly opening a 
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gap between it and the other states in the same shell. For example, in 234 Cm, this 
gap is more than 3.8 MeV, a value more than double of that observed in 218 Po. Thus 
in Cm, N n = 138 may act a subshell closure and the number of valence nucleons 
come out as N = 138. This cannot however explain why for Cf isotopes, it is 
necessary to have the number of valence neutrons at N n = N — 152. This may 
possibly be another effect of deformation. 

We expect that the observed smooth behaviour of the spectroscopic factor as a 
function of the parameter P should be reflected in some other properties also if the 
parameter P basically is a measure of proton-neutron interaction. With this in mind, 
we have chosen to study the ground state binding energy of the nuclei in this region. 
Our method is as follows. It is well known that correlations beyond mean field 
results are due principally to residual two body interaction. The residual interaction 
between similar nucleons is described by the zero range pairing force in the present 
calculation. However, the residual n — p interaction has not been considered in the 
present calculation. So for a chain of isotopes, the difference between the calculated 
and the experimental binding energies may be a measure of the strength of n — p 
interaction in a particular nucleus. However, not all differences can be ascribed to 
this effect. We have chosen one nucleus for each Z where, our scheme suggests that 
the number of valence neutrons is zero i.e. there is no valence n — p pair. We expect 
the effect of n — p interaction to be small in these nuclei and the difference between 
the calculated and experimental binding energies in these nuclei to be due to all 
the other effects combined. The difference in the change in the binding energy from 
the isotope with N n = for a particular Z between theory and experiment is taken 
as a measure of the contribution of N p N n interaction and expressed as A U7T . This 
difference has been plotted as function of P in Figure [3] in a scatter diagram. One 
can see a similar trend as observed in the case of spectroscopic factors particularly 
in the horizontal trend for low P values and the sharply upsloping values beyond 
P = 3. The results for Pu are somewhat ambiguous for two reasons. Firstly, as has 
already been pointed out, the N p and N n values are somewhat arbitrary. Secondly, 
the experimental binding energy value for Pu with N n = i.e. N = 126 is not 
known and we had to use an extrapolated value. In this case also, Pu nuclei do 
not conform to the general trend. However, the Cm isotopes also do not follow 
it. We have checked that the B(E2) values for 2 + — > + also show a similar type 
of behaviour, indicating that the quadrupole deformation starts to increase beyond 
P = 3, consistent with the fact that the proton neutron interaction is the major 
cause of deformation. 

To summarize, microscopic superasymmetric fission model has been employed 
to calculate the lifetime of alpha decay in even-even nuclei with Z = 84 — 98 and 
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N = 128 — 152. The interaction between the alpha particle and the daughter 
nucleus has been calculated in the double folding model using a density dependent 
NN interaction. The densities have been obtained using the Relativistic Mean Field 
formalism with the Lagrangian density FSU Gold. The spectroscopic factors for 
the decays are shown to vary smoothly as a function of effective numbers of valence 
nucleons, N p and N n . However, for a useful description, it was often necessary either 
to change the usual spherical subshell definitions or to utlize deformed subshells in 
counting N p and N n . The proton level nhg/2 and the neutron level 1^11/2 occur just 
after the closed shells at Z = 82 and N = 126, respectively and act sometimes as 
subshells. Binding energy systematics also follow a similar smooth curve indicating 
the role of the residual n — p interaction in spectroscopic factor. In view of the 
observed smooth behaviour, in future it will be interesting to investigate the origin 
of the observed effective valence nucleon numbers in more detail, particularly around 
Z = 94 and iV = 152.. 

This work was carried out with financial assistance of the Board of Research in 
Nuclear Sciences, Department of Atomic Energy (Sanction No. 2005/37/7/BRNS). 
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Table 1: Logarithm of experimental half life values and spectroscopic factors (S) of 
alpha decay obtained in the present calculation. 
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Figure 1: Spectroscopic factors S as a function of P = N p N n /(N p + N n ) for even- 
even nuclei between Z=84 and 98. The N p and N n values for different isotopes are 
explained in the text. 
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Figure 2: Calculated single particle proton levels near the proton shell closure Z=82 
for nuclei with iV — Z = 44. 
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Figure 3: The quantity A un plotted as a function of P for even-even nuclei between 
Z=84 and 98. See text for details. 
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